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Electronic Raman scattering measures a polarization-dependent scattering intensity which can 
provide information about the location of nodes in the energy gap of an unconventional supercon- 
ductor as well as its overall symmetry. In this paper, we calculate the Raman intensity in the 
presence of disorder for several models of the iron pnictide superconducting state. We include, for 
completeness, d-wave and isotropic s± responses in addition to more realistic extended s± supercon- 
ducting gaps. The effect of disorder is modeled using a self-consistent T-matrix approximation, and 
is studied in the limits of isotropic and intraband-only scattering. We show how recent experiments 
on Ba(Fei_ 2 ,Co x )2As2 may be consistent with "node lifting" by intraband disorder. 

PACS numbers: 



I. INTRODUCTION 



Electronic Raman scattering is sensitive to low-lying 
excitations and can be performed for various polar- 
izations of the incoming and outgoing photons. This 
polarization dependence preferentially samples different 
parts of the Brillouin zone, making Raman scattering 
an important tool with which to clarify the location 
of nodes and the symmetry of the gap in the super- 
conducting stated. Muschler et ali^ have recently per- 
formed Raman scattering measurements on single crys- 
tals of Ba(Fei_ a: Co 2 ;)2As2 (122) for two different con- 
centrations of Co. The data demonstrate the existence 
of low energy quasiparticle excitations in the supercon- 
ducting state. These quasiparticlcs could be a result of 
pairbreaking effects or due to the presence of nodes in 
the superconducting gap. In this paper we study the 
polarization dependent electronic Raman response in the 
presence of disorder in order to distinguish between these 
two possibilities. 

Thus far, other experimental probes have not pre- 
sented a convincing determination of a universal gap 
structure in the Fe-pnictides^. This may be due to 
variations in sample quality and resolutions issues of dif- 
ferent measurement techniques, but there is increasing 
speculation that these systems may possess a strong sen- 
sitivity of both the electronic structure and the pair state 
to small changes which affect electronic structure^. This 
is because the Fermi surface consists of several nearly 
compensated electron and hole pockets, and because the 
pairing state is probably of an extended- s type which may 
possess "accidental" nodes or deep minima, i.e. struc- 
tures depending on the details of the pairing interaction 
rather than the symmetry class. Nuclear magnetic reso- 
nance (NMR) studies^— showed a T 3 spin lattice relax- 
ation rate reminiscent of a gap with nodes. ARPES mea- 
surements on single crystals of 122-type materials^— 
measured the gap reporting isotropic or nearly isotropic 



gaps on all Fermi surface sheets. Penetration depth 
measurements^— have been fit both to an activated T- 
dependence, indicative of a fully gapped state, and low-T 
power laws, indicative of nodes in the superconducting 
gap. It is possible that these differences reflect genuinely 
different ground states in different materials due to in- 
trinsic differences in the pairing, but it is important to 
disentangle these effects from those which arise from dis- 
order, and in particular to distinguish between a disor- 
dered fully gapped state and a nodal state. 

A popular candidate for the ground state of ferrop- 
nictides is the so-called s± state proposed by Mazin et 
ali22. In this state, found within a simple spin fluctuation 
model with strong intcrband scattering between nearly 
nested electron and hole pockets, the gap is isotropic on 
both electron and hole pockets but with a sign change 
between the two. Further theoretical work^ir— consid- 
ered spin fluctuation pairing using band structures de- 
rived from tight binding fits to density functional theory 
results, and found s±-type states with highly anisotropic 
gaps, particularly on the electron sheets. We refer to 
this general class of states, which are in the same sym- 
metry class as ordinary s-wave and isotropic s± states, 
as "extended-s" states here. These calculations also re- 
ported the near degeneracy of d-wave gaps in certain sit- 
uations. The effect of nonmagnetic disorder on an s±, 
d-wave, and nodal extended-s is different and thus may 
help narrow the field of candidate superconducting gaps. 
In particular, it has been shown that for the accidental 
nodes of an extended s-wave state, intraband disorder 
can 'lift' the nodes, resulting in a fully- gapped stated. 
On the other hand, if disorder is of the intcrband type, 
low-energy impurity midgap states can be created^— 
similar to Yu-Shiba states due to magnetic impurities in 
ordinary superconductors. 

One complication is that changing the doping does not 
have the simple effect of creating more point-like disor- 
der. To what extent a dopant is charged or magnetic, 
long or short ranged, and what changes it makes in band 
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structure upon doping requires careful stud y 36 ' 37 and is 
difficult to include consistently in a theory capable of 
calculating observable quantities. In the most naive ap- 
proach, we will consider doping as loosely related to the 
concentration of scatterers in the sense that higher dop- 
ing is a dirtier sample, with the caveat that effects on 
the electronic structure and pair interaction may also be 
present. Here we report the results of a calculation of the 
Raman intensity for s±, c?-wave, and experimentally in- 
spired extended s-wave superconducting gaps, including 
disorder using a self-consistent T-matrix approximation 
(SCTMA). We show that for unitary scatterers, intra- 
band scattering will average the gap and 'lift' the nodes 
of an extended s-wave state, leaving a fully developed gap 
for higher scattering rates due to disorder. When a strong 
interband scattering component is present, the creation 
of an impurity band competes with the tendency to lift 
the accidental gap nodes. In contrast to the extended-s 
scenario, in an isotropic s± model, disorder causes the 
creation of low energy quasiparticles mimicking power 
law behavior of nodal states for some probes. One im- 
portant difference between such a sign-changing isotropic 
state and a nodal extended-s state, then, is that the ad- 
dition of disorder fills in the gap as opposed to possibly 
creating it. 

The form of the Raman spectrum will be influenced 
additionally by inelastic scattering processes. We expect 
these to be largely frozen out in the superconducting 
state at low energies, but they will be important for a 
proper treatment of the normal state spectrum and for 
energies near the maximum gap. Effects of this type will 
be treated elsewhere^ 

Our paper is organized as follows: in the first section 
we describe the theoretical background necessary to un- 
dertake this study. In the next section, we describe sim- 
ple one and two band results for the d-wave and s± cases 
to help ground our understanding. Finally, we show the 
results from studying anisotropic s-wave gaps on all four 
Fermi sheets which we suggest captures all the essential 
qualitative features of the real Co-doped I22s. 



II. THEORY OF ELECTRONIC RAMAN 
SCATTERING WITH DISORDER 

Raman scattering is the inelastic scattering of polar- 
ized light from a material (for a review see Ref(H). The 
cross section of the scattered light is proportional to^ 
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Here x 7 . 7 is the Raman effective density-density corre- 
lation function for symmetry channel 7. The vertex, 7k, 
accounts for the interaction of polarized light with charge 
density. The expression for the Raman effective density 
is P = X)k TkC^ o-Cicer- The full matrix Green's function 
in the presence of scattering in the superconducting state 
is 



UT + 6 k T 3 + A k T! 



e\ 2 



Ak 2 



where <D = to — £o> e k = £ k + S3, Ak = Ak + Si, and the 
Eq, are the components of the disorder self-energy propor- 
tional to the Pauli matrices r a in particle-hole (Nambu) 
space. 

Generally, the vertex is determined by both density 
and current matrix elements between the conduction 
band and the excited states. However, in situations 
where one is interested in qualitative results like the 
present one, model Raman vertices classified by symme- 
try can be employed. For a crystal with D^h tetragonal 
symmetry, in-plane charge fluctuations transform accord- 
ing to the full symmetry of the lattice (irreducible repre- 
sentation Ai g ), and do not change sign upon 90 degree 
rotation, or lower symmetry, as opposed to the B\ g and 
B 2g symmetry classes which change sign. Raman scat- 
tering probes long wavelength charge fluctuations. The 
Big and B 2g charge densities must average to zero within 
each unit cell, and are not coupled via the long-range 
Coulomb interaction. Thus the Bi g and B 2g Raman re- 
sponses for a multi-band system consist of the sum of the 
contributions from each band. However, A\ g fluctuations 
need not vanish over the unit cell, and therefore they can 
couple to isotropic charge density, giving rise the finite 
backflow. Therefore vertex corrections of both Coulomb 
and impurity type must be included in the A\ g channel, 
which therefore becomes extremely difficult to treat con- 
sistently with disorder. We do not treat this polarization 
channel in this work. Finally, we neglect direct inter- 
band and resonant contributions to the Raman vertices. 
In this approximation, the vertex for Raman scattering 
measures the effective mass around the FS. We use the 
notation that e is a polarization vector. 
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In a square lattice the B\ g and B 2g polarizations, the 
vertices can be expanded as: 



Bi g : 7k = cos(k x ) - cos(fc y ) 



where 



and 



X 7 . 7 M = / dre lUmT (T t/ 5 7 (t)/5 7 (0)) \ iUm ^ w+iS , (1) 
Jo 

and 

X (q -> 0) = r^^Tr( 7 ^r 3 G(k,ia;„)T3G , (k,ia; n +ffl ro )) 

n k 
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7k = sin(fe K ) sm(k y ) 



which along a circular Fermi surface become cos(20) and 
sin(20), respectively. It should be noted here that in the 
iron-pnictides use of the 1-Fe vs. the 2-Fe Brillouin zone 
will interchange B\ g and B 2g symmetries. 
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To include the effect of disorder we sclf-consistently 
solve for the self energies by including all scatterings off 
a single impurity and performing a disorder average. This 
is represented diagrammatically in Fig. [TJ The T-matrix 
can be defined as 



G — Gq + GqTGq + 



(2) 



+ 



+ 



+ 



FIG. 1: Diagrams representing the self-consistent T-matrix 
approximation (SCTMA). The dashed lines are scattering off 
an impurity, and the single line is the self-consistent green's 
function. 

Using the Nambu notation, the single-particle self en- 
ergy in a superconductor can be decomposed as: 



S(k,w) 



E a (k,w)f Q , 



(3) 



where f a are the Pauli matrices and fo is the unit ma- 
trix. Note that the band index is implicitly contained in 
the k index since we restrict pairing to individual Fermi 
surface sheets. Treating impurity-scattering in T-matrix 
approximation gives rise to the following self-energy 



E(k,w) = niTkk(w), 



(4) 



where n., is the impurity concentration and Tkk{w) is the 
diagonal element of the T-matrix 



Ikk'M = Vwh + VW'Wk", w)2W(w) • (5) 

k" 

We define 
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where Ui is the density of impurities, n of electrons, and 
Nq the density of states at the fermi level. For a constant 
potential, the case which we will consider in this paper, 
this expression becomes the series: 

E(w) = mVofs G ( k '' u ) V o^T 

n k' 

Later we will restrict Vkk" to be constant for particu- 
lar sets of momenta, either to allow transitions between 
all the Fermi sheets or to restrict transitions to remain 
within Fermi sheets. 

The self-energy E(fc,w) has to be solved self- 
consistcntly in combination with the single-particle 
Green's function 



After solving for the self-energies, we insert them into 
the general expression for the Raman response. Begin- 
ning with a spectral representation: 



1\ 7mG(k,x) 



G(k, ioj n ) = / dx ( — 

J \ 7T J lW n — X 

we arrive at a zero-temperature long-wavelength form, 
f° 1 

Imx-y Jil) = V / dx-jlTr[ImG(]i, x)r 3 ImG(k, x+n)r 3 ] 
k J n 

In terms of retarded and advanced Green's functions, 



Jmx 7l7 (n) = ^(N(<f>)Im I d^l I 



[f hh {x - - f (x - n, x ) 



F AH (x -n,x)+ F [x - n, a;)]) 



where 



G^w)" 1 = G (k,( 



E(w). 



(6) 



F a ' b = Tr[G a (k, x - ft)r 3 G b (k, x)t 3 ] a,b = A,R 

We have taken N(<f>) = Nq in this paper. The angular 
brackets denote an average over the angle <f> around the 
Fermi surfaces. 

We note here that for crossed polarizations selected 
Big and B 2g channels, the Raman response is addi- 
tive for a many-sheeted Fermi surface. We also note 
that in the present approach we have neglected T-matrix 
impurity vertex corrections for the Raman response. 
While for general momentum-dependent Ti^k', impu- 
rity vertex corrections are necessary for all channels, 
for momentum-independent scattering vertex corrections 
do not essentially modify the response for B\ g and B 2g 
polarizations^. 



III. RAMAN SCATTERING FOR DIRTY 
d-WAVE AND s±-WAVE CASES 

First we present results for d-wave and s±-wave super- 
conducting gaps. Early on, the s± was proposed as a 
candidate for the gap structure in the iron-pnictides^. 
Indeed, for certain materials parameters, in particular 
the case where a Fermi surface pocket at (7r,7r) appears, 
this state is crudely consistent with multiorbital spin- 
fluctuation calculations^^—. The qualitative features 
of the Raman response for such a state can be under- 
stood from a minimal two-Fermi-sheet model. 

NMR and penetration depth measurements in some 
Fe-pnictidc materials displayed power-law temperature 
dependences suggestive of nodes in a single band 
model£~— ii£r— . The phenomenology was similar to what 
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ers, a low-energy impurity band is created in the density 
of states, show in Fig. O The clean Raman response 
simply reflects two clean s-wave gaps with a sharp gap 
edge at 2Ao- Note the line-shape in a fully gapped super- 
conductor. This qualitative feature seems differs from a 
nodal supcrconductors's line-shape, which still possesses 
a peak, only more symmetric about the center energy at 
2Ao- One can understand the features in the Raman re- 
sponse as crudely similar to a convolution of the density 
of states with itself, so the effect of the impurity band on 
the Raman response is to create a nonzero threshold at 
uj < 2Aq, as shown in Fig. [3J corresponding to transi- 
tions between the impurity band and the gap edge. There 
is, in addition, a very small contribution from scattering 
within the impurity band itself at low energies (see insert 

Fig. EH. 



FIG. 2: (Color online) Quasiparticle density of states N(u>) 
for isotropic s± state, normalized to normal state density of 
states A^o vs. u/Ao, where ±Ao is the value of the gap on 
hole and electron sheets. Nq is assumed constant on all Fermi 
sheets. Shown are various interband impurity scattering rates 
r in units of Ao. 
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FIG. 3: (Color online) Raman response Im x( w ) vs - w/A for 
both Big and £>2g polarizations for an s± state as in Fig. [2] 
for various interband impurity scattering rates T, in units of 
A . 

was seen in cuprates, so it is useful to include a study of 
the d-wave case for comparison, to how the presence of 
a node in the superconducting gap manifests itself in the 
Raman response. The d-wave Raman response has been 
previously studied, for example in Refs. li^liil andl43l 

We first consider the s± state on two circular Fermi 
sheets, each with an isotropic constant gap in the su- 
perconducting state A± = ±Arj. The gaps differ by a 
minus sign, so inter-band scattering gives rise to pair- 
breaking and violates Anderson's theorem. As we in- 
crease the scattering rate due to strong isotropic scatter- 
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FIG. 4: (Color online) Density of states N(u)/Ao vs. en- 
ergy oj/Ao for a d-wave superconductor for various values of 
scattering rate T/Ao in unitarity limit. 

In Figs. 0] and [SJ we show the corresponding quanti- 
ties calculated in the simple 1-band d-wave case, A(6) = 
Aocos(2#), as a canonical example of what one expects 
for a nodal unconventional superconductor. One funda- 
mental difference with the s± case is the presence of low 
energy quasiparticles in the clean limit. The nodes in 
the superconducting gap allow for excitations in the low- 
energy density of states, shown in Fig. UJ Furthermore, 
the Raman polarization which weights the nodes (i?2g, 
in the d-wave case) yields differing power laws in uj from 
the polarization which vanishes along the nodal direc- 
tions (Bi g ), in contrast to the isotropic case where both 
of the Raman responses (B\ g and i?2g) are the same. In 
the d-wave case, this is signaled by the presence (Bi g ) or 
absence (i?2g) of a large peak at 2Ao- In both cases, in- 
creasing the scattering rate due to disorder increases the 
size of the impurity band, but the qualitative features of 
the d-wave Raman response, shown in Fig. [5] (a peak at 
2Ao, excitations down to uj = 0, differing response for 
differing polarizations) are unaltered by disorder. Be- 
yond blurring of the sharp features, the effect of disorder 
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FIG. 5: (Color online) Effect of disorder on T = Raman 
response of d wave state vs. energy cj/Ao. Shown are two 
polarizations, Bi g and -B29, for various values of scattering 
rate T/Aq in unitarity limit. 
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FIG. 6: Fermi surface from a five-orbital tight binding model 
with parameters chosen to match the density functional the- 
ory result. 



will change the low frequency behavior between from w 3 
to lu in the B\ g polarization for the ci-wave case 3 -^. 



IV. RAMAN RESPONSE IN MODELS OF 
FERROPNICTIDES 



The full band structure of the Fe-pnictide LaOFeAs, 
determined by density functional calculation s 44 ! 45 , can be 
accurately parameterized by a tight binding model with 
S-bands^ 4 -. Generally, these materials have four Fermi 
surface sheets, shown in Fig. |6] in the "unfolded" or 
1-Fe zone: two hole pockets around the gamma point 
and two electron pockets. The two Fermi surfaces about 
the gamma point are referred to as a sheets and the 
other two are /? sheets. The exact details of the band 
structure are sensitive to doping 5 -, and an additional hole 
pocket around the (n, tt) point can occur. We do not 
include the (it, 7r)-pocket because its effect is expected 
to stabilize isotropic gaps in the context of multiorbital 
spin-fluctuation calculation o 5 ' 27 i 29 , and because it occurs 
(within a rigid band shift implementation of doping) for 
the hole-doped cases only. We focus here on electron 
doped materials on which experiments have been per- 
formed, which appear to show nodes or deep gap minima. 

To model the Ai g , B\ g , and B 2g polarizations mea- 
sured in Muschler et alj^ the same choices were made 
for the extended s-wave gap on the a sheets (around the 
r point) and /3 sheets (around the M point) as in that 
paper: 
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(7) 
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The vertices for this state arc chosen to be: 



e*i 7i g = 0. 7 2g = 0. (11) 

a 2 -fi g = .25(-2) sin(t9) cos(t?) l2g = .25 cos(20)(12) 

ft 7i ff = .5(-2)sin(0)cos(0) l2g - +1. (13) 

P2 yig = .5(-2)sin(0)cos(0) 72g = -1. (14) 

Note that the notation is now in the 2-Fe zone, for easy 
comparison with Ref. The vertices have been chosen 
to model experimental results with several constraints 
in mind. First, we must respect the underlying symme- 
try of the polarization state. In the B\ g polarization, 
the response weighted is largely away from any Fermi 
sheet, which is reflected by a flat response in the model 
and data. The B 2g polarization samples the electron (/?) 
sheets^. In the data and the model, there is a strong 
T-dependence to the data and a peak which appears be- 
low T c , reflecting the corresponding gap function on these 
Fermi surfaces. More realistic calculations of the vertex 
functions 7k will be necessary for quantitative compari- 
son with experiment. 

Figure [7] shows the clean results for this model. 
In Muschler et alj^, single crystals of Co doped 
Ba(Fe 1 _ 2: Co :I ;)2As2 with x = 0.061 and 0.085 were stud- 
ied. First, at x=0.061, the B 2g spectra in the supercon- 
ducting state are strikingly different from the B\ g and 
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A\ g spectra, which is not possible for an isotropic gap. 
There is a sharp peak in the B 2g polarization whose shape 
is not asymmetric as would be characteristic of a full gap, 
already suggesting the presence of nodes. The B\ g polar- 
ization shows almost no change upon entering the super- 
conducting state, which is consistent with vertices which 
probe the regions of the Brillouin zone without a Fermi 
surface. For Ba(Fci_ T Co x )2As2 with x=0.061, there is a 
nonzero Raman intensity down to zero frequency, which 
indicates the presence of low energy quasiparticles. Thus 
far, this is reminiscent of the one-band d-wave model. 
To explore this further the low frequency behavior of the 
Raman intensity can be examined. It is straightforward 
to show that the observed w 1 / 2 -dependence of the Bi g 
is characteristic of a marginal or "kissing" node, where 
the node just touches the Fermi surface. Examination at 
other dopings can now shed light on the nature of the 
state. A revealing aspect of the data£ is the x=0.085 Bi g 
spectra in the superconducting state. For this higher 
doping a finite gap of 10 cm -1 can be resolved. The 
changing structure observed at two dopings suggests that 
impurities can play an important role. 



Clean B2g 

Clean B1g 




FIG. 7: (Color online) Clean Bi g and B29 spectra for the 
model of the superconducting state defined in Eq. [7lll0£ and 
in Eq. ©-COB • 

The effect of disorder was calculated using the SCTMA 
in two limits: isotropic intcrband scattering for which 
V(k, k') = Vb, and intraband-only scattering for which 
V(k, k') = if k, k' are on different Fermi sheets and Vb 
if k, k' are on the same Fermi sheet. It was assumed that 
the potential Vb is large (unitarity limit), simulating the 
large density functional theory effective potential found 
for Co dopants in Ba-123^. In Fig. [8] we can see the 
effect of increasing the intraband scattering rate on the 
density of states. The gap A^^f?) has been chosen so 
that in the clean limit it has a marginal or "kissing" 
nodes leading to a characteristic lo 1 ^ 2 behavior in the 
Raman respons e 2 ' 41 , and the asymptotic low-w behavior 
of N(ui) is indeed also a; 1 / 2 , as seen in Fig. [5J however 
there arc additional features at low energies due to the 



small minimum gaps on the hole sheets. These features 
are largely suppressed in the clean Raman response, as 
shown in Fig. [7j In the presence of intraband disorder, 
the nodes on the (3 sheets are lifted immediately, as seen 
in the density of states, creating a small but complete 
spectral gap A m i n , as discussed in Mishra et al<^ This, 
in turn, is reflected with the creation of a gap in the 
Raman intensity in Fig. [5J where there is a disorder 
induced gap edge in the response up to w = 2A m i n . 



2.5 




0.5 1 1.5 2 2.5 3 3.5 

<o/a 

FIG. 8: (Color online) Density of states N(uj)/N for the 
extended s-wave state given in Eqs. Eq. I7I10I vs. u /An for 
unitary intraband scatterering rates r/An. Insert shows low- 
energy behavior. 
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FIG. 9: (Color online) Raman intensity for the extended s- 
wave state given in Eqs. Eq. I7I10I vs. a;/ An for various 
unitary intraband scatterering rates r/An and polarization 
states Bi s and £?2 S in the 2-Fe zone. Insert: low energy region. 

The isotropic scatterers differ from the intra-band scat- 
terers fundamentally in that inter- Fermi-sheet transitions 
allow for significant pair breaking, since the extended s- 
wave state changes average sign between electron and 
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hole Fermi surfaces. So in addition to the gap-averaging 
that occurs as we scatter from k to k', there is the addi- 
tional effect of creating a low energy quasiparticlcs "im- 
purity band." This effect prevents the lifting of nodes, 
so while the gap becomes more isotropic, the impurity 
band is the most important low energy effect in the den- 
sity of states. Indeed, as shown in Fig. [TU1 the DOS 
is reminiscent of the dirty d-wave case. There is also a 
striking similarity in the Raman intensity shown in Fig. 
QT]and the rf-wave case. There are low energy quasipar- 
ticles all the way down to zero frequency, and we observe 
an enhancement of low frequency spectral weight. 
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FIG. 10: (Color online) The density of states for an extended 
s-wave state given in Eqs. Eq. 171101 vs. w/Ao for unitary 
isotropic scatterering rates F/Aq. 
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FIG. 11: (Color online) Raman intensity for an extended s- 
wave state of Fig. [TO]for various unitary isotropic scatterering 
rates T/Ao and polarization states Bi g and B2 g in the 2-Fe 
zone. Insert: low-energy region. 



CONCLUSIONS 



Measurement of the Raman spectra in 
Ba(Fci_ :c Co :c )2As2 revealed that at roughly 8% doping 
there is a gap in the response up to approximately 10 
cm , in contrast to the 6% doped sample where there 
is a nonzero response down to zero frequency. This 
intriguing feature suggests a relation between cobalt 
doping, disorder, and the structure of the underlying 
superconducting gap. We have presented calculations of 
the Raman response as a function of frequency for vary- 
ing scattering rates due to disorder in the self-consistent 
T matrix approximation. First, we demonstrated the 
differing effects of disorder on s± and d-wave gaps to 
illustrate how a nodal superconductor contrasts with 
a sign-changing gapped state. The absence of low 
energy quasiparticlcs is an indication of a fully gapped 
state, but for systems with interband disorder scattering 
an impurity band is created which can mimic nodal 
behavior. We found, however, that in the presence 
of disorder the form of the main 2A peak is largely 
unaffected; thus its observed symmetric form in the 
Muschlcr ct al experiment is suggestive evidence for true 
nodal or near-nodal behavior. 

More definitive information is provided by the Raman 
response at low energies. To study the anisotropic s-wave 
state in more detail, we have presented a set of model 
gap functions on the 2D Fermi surface sheets which cap- 
ture the essential features of the experimental Raman 
response^. It is immediately clear that for the dopings 
where measurements have been performed a nodal state 
is likely because of two factors: low frequency power law 
behavior of the Raman intensity down to zero energy, 
together with the very different responses for different 
polarizations. We then examined different kinds of dis- 
order on such a state. A remarkable effect is the lifting 
of nodes by intraband scatterers for a highly anisotropic 
s-wave state, in rough agreement with the qualitative 
effect of gapping the Raman spectrum observed as the 
doping increased from 0.061 to 0.085 in Muschler et al. 
in Ba(Fei_ x Co x )2As2 . The predicted effects of strong 
interband scattering at low energies were not observed, 
implying that the Co in this system acts primarily as an 
intraband scatterer— . Of course, the present calculations 
include only the effect of the impurity on pairbreaking, 
and neglect any direct effect on the pairing interaction, 
so direct conclusions may not be drawn. 

This feature distinguishes nodal anisotropic s-wave 
states from isotropic s± states in terms of how the 
response, and underlying superconducting gap, evolves 
with increasing disorder. It also is not possible in non- 
Ai g representations of the gap, like d-wave, because the 
nodes cannot be removed for lower symmetry supercon- 
ducting gaps. We believe that these distinctions will 
provide a useful basis for the interpretation of Raman 
scattering experiments as doping or irradiation is tuned 
systematically in different samples. 
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